A new model of turbulent relative dispersion: a self-similar telegraph equation based 

on persistently separating motions 
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Turbulent relative dispersion is studied theoretically with a focus on the evolution of probability 
distribution of the relative separation of two passive particles. A finite separation speed and a finite 
correlation of relative velocity, which are crucial for real turbulence, are implemented to a master 
equation by multiple-scale consideration. A telegraph equation with scale-dependent coefficients 
is derived in the continuous limit. Unlike the conventional case, the telegraph equation has a 
similarity solution bounded by the maximum separation. The evolution is characterized by two 
parameters: the strength of persistency of separating motions and the coefficient of the drift term. 
These parameters are connected to Richardson's constant and, thus, expected to be universal. The 
relationship between the drift term and coherent structures is discussed for two 2-D turbulences. 

PACS numbers: 47.27.Qb, 47.27.-i, 05.40.-a, 92.10.Lq 



Turbulent transport and mixing underlie wide range 
of phenomena from star formations to coffee in a cup. 
However, the mechanisms of their significant deviation 
from molecular counterparts are not well understood yet 
even in such a simple case as relative dispersion of passive 
particles. 

In the inertial subrange, reflecting the scaling law of 
turbulent velocity fluctuation, well-known Richardson's 
t 3 law is realized 01 • F° r the probability density func- 
tion (PDF) of separation r, P(r,t), quite a few mod- 
els and theories consistent with Richardson's law have 
been proposed 0]. Besides recent progresses of particle 
tracking techniques and numerical simulations enable di- 
rect experimental investigations of the separation PDF 
[EL IE • Although accuracies of the experiments arc not 
enough, most of their results are so far close to the pre- 
diction of Richardson's diffusion equation 0: 
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where K(r) cx r 4 / 3 in the Kolmogorov scaling |3( and d 
is the spatial dimension. This closeness indicates Eq. (Tj^) 
is a good basis for description of turbulent relative dis- 
persion. Eq. (Tj^) ought to be exact only if the relative 
velocity is (5-correlated in time 0,0- However, there are 
spatio-temporal correlations in turbulent flows as implied 
by existence of coherent structures. This indicates rela- 
tive velocity cannot be (5-correlated in time. To resolve 
the inconsistency, Eq. (Tj^) has to be extended to include 
these correlations. As mentioned in detail below, we treat 
these correlations as those not in time but in scale-space. 
That is, we focus on correlations between scales r and pr, 
where p is a scale multiplier. This treatment is appro- 
priate for describing self-similarity. Employing multiple- 
scale consideration with correlation in scale, we derive a 
telegraph equation with scale-dependent coefficients, Eq. 
(|f and obtain a similarity solution, which never co- 
incides with Richardson's one markedly in the tail part 



even in the long time limit. 

Coherent structures observed in turbulence must share 
their origin with finite correlation and self-similarity, so 
that particles appear to be advected according to coher- 
ent structures. For example, in two-dimensional (2D) in- 
verse cascade (IC) turbulence, particles separate step-by- 
step through nested cat's eye vortices with scattering by 
stagnation points ! m 2D free convection (FC) tur- 
bulence, particles separate through advection by stretch- 
ing and folding plumes (Fig. Q j| . Clearly, separation 
processes in these two systems are different. However, 
despite these differences, effects of coherent structures 
on relative dispersion appear in the same way, i.e., per- 
sistent separation: persistent expansion and compression 
of a relative separation. Sokolov et al. model such mo- 
tions based on Levy walk. Their model consists of persis- 
tent separation ceased by probabilistic turn in direction 
0, 0| . They introduced also the persistent parameter 
P s , the ratio of the correlation length to the scale. 

The telegraph model was introduced to implement (I) a 
finite diffusion speed, and (II) a finite correlation time to 




FIG. 1: A snapshot of vorticity field of 2D-FC turbulence. 
Shading represents intensity of vorticity. The scale denoted by 
the straight line is WOr/g and in the inertial subrange, where 
r/e is the thermal Kolmogorov scale 0]. 
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the diffusion process 0, 0] . It has been widely applied 
to various diffusion phenomena from molecular diffusion 
[l3| to population dynamics (lif. 

To satisfy the self-similarity of turbulence, (I) and 
(II) are extended to be scale-dependent according to 
the following scaling assumptions: v(r) = Ar 1 ^ 9 and 
T c (r) = r/v(r) = A~ x r 9 , where v(r) and T c (r) are rela- 
tive velocity and characteristic time, A a constant [lij . 
and g a scaling exponent . 

First, we focus on a scale f in the inertial subrange. 
The correlation length and time for a separation to be 
expanded or compressed persistently are defined as P^f 
and PgT c (f), respectively. We assume Pf- is the order of 
unity, and then, p out = 1 + Pp. We call this scale outer. 

We consider the evolution of probability density of the 
separation in a small region around f where its extent is 
much smaller than P^f. Since P^ f and P±T c (r) are re- 
garded as constants, we can apply the approach deriving 
the telegraph equation to this small region We call 
this scale inner. 

We divide the inner region into shells defined as 
[r n ,r„+i) where r„ = f£™, n is integer, and £ is close 
to unity, that is, p m = £™. Since pi n — 1 <SC p ou t — 1, we 
deal with two scales in scale-space. To take into a finite 
speed, pass-through time r„ of the n-th shell, the time 
for a relative separation to expand (compress) through 
the n-th shell, is defined as follows: 



Ar» 
v(r n ) 



v(r n ) v(r n y 



(2) 



where 7 = log£ (7 <C 1), and v(r n ) is the relative velocity 
at a spatial scale r n . This relation means 0(t„) = 0(7), 
which is the key difference from the diffusion equation 
case including Richardson's equation. 

We introduce probabilities Q^: the probability for a 
relative separation to be expanding (Q+) or compressing 
(Qn) m the n-th shell. Transition probabilities Ap (r n ), 
the probability from expansion to compression (+) and 
the opposite (— ) during r„, must be self-similar in outer 
scale. Therefore the simplest form of Ap ± (r„) is 



Ap ± (r„) = 



PpT c (f) T c (f) 



(3) 



where P s ± T c (f) are correlation times depending on the 
direction and X ± = 1/P±. Note that Pf~T c (r) are outer 
scale and, thus, constants in the inner region. This form 
was first given by Sokolov et al. 0] • Unlike their model, 
we assign different values to X^ 1 to represent difference of 
persistency between expansion and compression. 

Based on above considerations, we construct master 
equations for Q„~(t) and Qn(t) as follows: 
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The last terms in the r.h.s. of Eqs. 0} and JSJ denote the 
remainder of particle pairs leaving the n-th shell during 
t„_i and t„ + i, respectively. Because Ar„ <C 1 and r„ <C 
1, (i) the Kramers-Moyal expansion and (ii) addition and 
subtraction of Eqs. 0} and © lead to 
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In order to combine Eqs. JB} and Q, first, divided by r n 
and differentiated with respect to n, Q leads to 
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where Q„ = Q+ + Q~, A = A+ + A~, and S = X + — A - . 
Next, differentiating © with respect to t and substitut- 
ing (JH} to it, we obtain an equation for Q n around f: 



T c (f) d 2 Q n , aQ r . 
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Taking the limit as 7 — > 0, r n approaches r. 
And thus we obtain the following relations: Q n (t) = 
fr" +1 P(r,t)dr w jrP(f,t), d/dn w jrd/dr, and 7/t„ « 
v(f)/f. Replacing f with r, we finally have an equation 
for P(r,t): 
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Then Eq. (|10|l with the scaling assumptions yields 
T c (r) d 2 P i dP 

a 
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a—[v(r)P], (11) 



where D (r) is Richardson's diffusion coefficient, 
AA -i r 2- fl) and a = ^ -2g + (S)A- 1 = ctA- 1 . The 
conservation of probability is satisfied at least for the 
similarity solution. 
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In the limit of infinite speed and ^-correlation, the first 
term of Eq. (|If I) disappears, which is the same form as 
Palm's equation (see p. 575 of [2(J- Non-Richardson 
terms, the first term in the l.h.s. and the last one in the 
r.h.s. of Eq. (fill) , describe effects of persistent separation. 
The last term in the r.h.s. of Eq. I|llfl is a drift term con- 
sistent with the scaling assumptions. The drift velocity 
is —ov{r) and, hence, the direction is determined by a 
which consists of the "scaling-determined" part, d — 2g, 
and the "dynamics-determined" part, S. 

Now, let us calculate the similarity solution of Eq. (|f f [) . 
To obtain it, we introduce the similarity variable £, 



9 \ 
— r~ 9 r 9 
(r 2 )V2/ °« At * 



(12) 



where we set (r 2 } 1/2 = C R {At/X)^ 9 21]. Substituting 
P(r, t) — F(^)(r 2 } -1 / 2 into Eq. fTT|) . we obtain an equa- 
tion of F(0 = F{ri/C 9 R ) = F(rf): 



drj 2 

where r\ — C'^ 9 ^, 
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and d = d — a . Eq. (|13fl has three fixed singular points: 
£ = and £ = £± = ±Xg/C R . Starting from an initial 
condition, where all particles occupy a small and compact 
region, the solution is confined in [0,£+]. This finitcness 
of the similarity solution is due to the assumption of a 
finite speed and, thus, the maximum relative separation, 
W - e+{r 2 ) 1/2 = (Agt) 1 / 9 , exists 0. 

To see asymptotic behaviors of the similarity solution, 
we consider two limits: the small (£ <C 1) and the max- 
imum (£ ~ £ + ) separation regimes. In the small sep- 
aration regime, £ -C 1 or r <C (r 2 ) 1 / 2 , we obtain the 
following solution: 



F(0 oc {MY 2 exp(/3i0, 



(16) 



where (3\ = (d — o — l)/g and 02 = C 9 R / g 2 . This solution 
is also a similarity solution of Palm's equation, because 
the first term in Eq. (|llfl can be neglected under the 
condition r <C (r 2 ) 1 / 2 . We call this limit the diffusive 
regime. 

On the other hand, in the maximum separation regime, 
the frontal edge of the PDF is abrupt at r max . The func- 
tional form of the edge is approximated to the first order 
by 



TABLE I: Estimated values of 8 and a in 2D-IC and -FC 
turbulences. 



Type 


Method 
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a 


2D-IC a 


Exit-time PDF 


-1.48 


-0.81 


2D-IC i ' 


Separation PDF 


-0.87 


-0.20 


2D-FC C 


Exit-time PDF 


-0.71 


0.49 



"From DNS results by Boffetta and Sokolov 
''From DNS results by Goto and Vassilicos 
c From our DNS with resolution 2048 2 [Trt . 
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We call this limit the telegraphic regime. From this 
expression, it is clear that if A < d, P{r, t) — > oo as 
r — ► r max . That is, most of particle pairs are accumu- 
lated at the frontal edge, where relative separations ex- 
pand away without changing their directions. However, 
this situation is unrealistic in real turbulence, so that A 
is considered to be greater than d. 

Control parameters of our model are A and a. The 
functional form of the PDF, F(£), is determined by these 
parameters: a controls mainly the diffusive regime, and 
A does the telegraphic regime. As A = 1/-P+ + 1/-P S ~, 
it represents the strength of persistency of moving direc- 
tion. On the other hand, as o\~ x is the coefficient of the 
drift term, a represents total effects of persistent sepa- 
rations and probabilistic transitions. Because persistent 
motions model the advection by coherent and self-similar 
flows, a seems to characterize the average effects of flow 
structures such as coherent structures on dispersion pro- 
cesses. In order to calculate the value of a, we have to 
estimate "dynamics-determined" part of it, S. 

To estimate S from direct numerical simulations 
(DNS), we use the PDF of exit-time @, El- Exit-time 
for particle-pairs experienced many turns form an expo- 
nential tail (for details see 0). In our case, the slope 
is evaluated with Eq. (|llfl in the limit of infinite time, 
i.e., Palm's equation, where the slope is related to 6 [23j |. 
We can also estimate S directly from the separation PDF 
around r < (r 2 ) 1 / 2 . In the 2D-IC case, Goto and Vas- 
silicos estimated a = 2/(g — 5) by fitting the similarity 
solution of Palm's equation to the PDF jjj. 

Table H] shows the estimated values of S and a for 2D- 
IC and -FC turbulences. These results indicate that the 
drift term of (jl 1|> enhances diffusion in the 2D-IC case but 
suppresses diffusion in the 2D-FC case; Compression of 
relative separations in 2D-IC turbulence but expansion of 
them in 2D-FC turbulence are comparatively restricted, 
respectively. This remarkable feature is considered to 
be induced by the difference in flow structures between 
2D-IC and -FC turbulences: "cat's eye in a cat's eye" 
structures 7] and string-like structures 0- We, there- 
fore, expect that a can characterize coherent structures. 

In Fig.El the similarity solution F(£) obtained numer- 
ically for various values of A is shown. A cut-off scale 
corresponding to r max can be seen. Our similarity solu- 
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FIG. 2: PDFs of relative separation for the 2D-IC case. Solid 
lines represent similarity solutions of our model with A = 4, 
6, 8, 10, and 15 from the curve with the smallest cut-off scale 
to that of the largest. The dashed line represents A = oo, or 
the similarity solution of Palm's equation 11161 . a is fixed for 
1/3: d — 2, g — 2/3, and 8 — — 1. The inset is the same plot 
in linear scale. Lines correspond to A = 4, 6, 8, 10, 15, and 
oo from the curve with the lowest peak value to that of the 
highest. 

tion approaches Palm's one as A gets larger. However, 
even for large A, the difference between them in the tail 
part is so evident that effects of persistent separation are 
not negligible. 

In summary, we derived a telegraph equation with 
scale-dependent coefficients, into which finite separation 
speed and self-similarity are incorporated, by employing 
multiple-scale consideration in scale-space. Then we ob- 
tained a similarity solution of it. In the diffusive regime, 
£ <C 1, the similarity solution coincides with that of 
Richardson's diffusion equation with the drift term, i.e., 
Palm's equation; in the telegraphic regime, £ ~ £ max , the 
finiteness of separation speed is realized and the separa- 
tion PDF is abrupt at r max . Therefore, finite separation 
speed is crucial for description of the tail part of the 
separation PDF unless relative velocity is ^-correlated in 
time. 

The drift term of Eq. is induced by the devia- 
tion of the difference of persistency between expansion 
and compression of relative separation, S, from "scaling- 
determined" value, 2g — d. The direction of the drift is 
determined by the sign of —a = (2g—d)—S. We estimated 
the value for two 2-D turbulences, inverse cascade (IC) 
and free convection (FC), and found that positive and 
negative drift is imposed in the 2D-IC and -FC cases, re- 
spectively. We conjecture that this remarkable difference 
corresponds to the different types of coherent structures 
in the background flow. We need more precise investiga- 
tions to obtain an evidence for this. 

We neglected two significant effects: distribution of 
separation speed and intermittency. However, intermit- 
tency of relative velocity is negligible in the two 2D tur- 
bulences dealt with in this letter. Besides, not the distri- 



bution but the finiteness and self-similarity of separation 
speed is crucial for the existence of the cut-off of the sep- 
aration PDF. We are expecting experiments which can 
resolve the tail part of the separation PDF. 
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